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Abstract 

In this paper, we perform a comprehensive study to probe the effects of large extra di¬ 
mensions through double Higgs production in proton-proton collisions at the center-of-mass 
energies of 14, 33, 100 TeV. We concentrate on the channel in which both Higgs bosons decay 
into bb pair and take into account the main background contributions through realistic Monte- 
Carlo simulations. In order to achieve an efficient event reconstruction and a good background 
rejection, jet substructure techniques are used to efficiently capture the boosted Higgs bosons 
in the final state. The expected limits on the model parameters are obtained based on the 
invariant mass and the angular properties of the final state objects. Depending on the num¬ 
ber of extra dimensions, bounds up to 6.1, 12.5, 28.1 TeV are set on the model parameter at 
proton-proton collisions with the center-of-mass energies of 14, 33, and 100 TeV, respectively. 
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1 Introduction 


The SM-like Higgs boson discovered by the CMS and ATLAS jT, 2] experiments at the LHC 
indicates a strong evidence for the proposed mechanism of the spontaneous electroweak symmetry 
breaking (EWSB) in the Standard Model (SM). However, further efforts are ongoing to test the 
characteristics of this newly observed particle against the SM Higgs [3] in terms of its properties 
and couplings to the SM particles. Another important question is the EWSB behaves like what 
is predicted by the SM whether or no. To answer this question the Higgs potential has to be 
examined up to higher order through measurements of Higgs boson self-couplings which has many 
interesting phenomenological implications ® u 0 El i. One way to explore the Higgs self¬ 
couplings is via measurements of the di-Higgs and triple-Higgs productions at the LHC and future 
planned hadron colliders [9j. 

On the other hand multiple Higgs boson production also could be used to search for new 
physics especially in the high invariant mass region either for resonant or non-resonant possible 
physics beyond the SM. In the resonant case, there are several interesting hypotheses which permit 
new resonance decaying to Higgs pair such as Randall-Sundrum radion m and CP-even heavy 
Higgs in the next to minimal supersymmetric standard model(NMSSM) |11| . In addition to that, 
there are already several studies to search for new physics through di-Higgs final state which can 

be found in na [la nn nu nsi nzi nn ns] . 

In the non-resonant searches, a possibility is to use the di-Higgs final state at the LHC and 
Future Circular Colliders (FCC) to search for large extra dimensions in Arkani-Hamed, Dimopou- 
los, Dvali (ADD) scenario [20]. They proposed the large extra dimension scenario as a solution 
for the hierarchy between the scale of electroweak and Planck scale [201 ED 122] . According to 
their model, the SM fields are confined to the 3+1 space-time dimensions while gravity can freely 
propagate into the multi-dimensional space-time 4 + tied, where tied is the possible number of 
extra dimensions. This leads to propagation of gravity field flux into the entire 4 +ued dimensions 
which leads to dilution of the power of gravity in the common 3+1 dimensions. The reduction 
of the gravitational flux can be quantified by applying the Gauss’s law. The result expresses the 
relation between the ordinary fundamental Planck scale Mpi in 3+1 common dimensions and 
Planck scale in 4 + tied dimensions denoted by Me according to the following relation: 


M ™ ED+2 


rsj 


Mpi 

Rn E D 


(1) 


where R is the size of extra dimensions. According to the ADD model motivation, if one assumes 
the Me ~ Mew ~ 1 TeV the size of extra dimensions for tied = 2 to 7 can be varied from few 
centimeters down to few femto-meters. 

Based on the ADD scenario, many phenomenological studies have been preformed to find the 
possible observations in the particle colliders mmm- In these works, graviton in the multi¬ 
dimensional representation equivalently interpreted as towers of massive Kaluza-Klein (KK) modes 
or Gfefcwhich can couple to the SM particles through SM energy-momentum tensor. The resulting 
effective model provides different experimental signatures such as virtual exchange of graviton and 
direct graviton emission at colliders. 

Although the coupling of each KK mode G^k with the SM gauge bosons or fermions is 
suppressed by the Planck scale, the summation over all KK-modes with tiny mass splitting 
Auikk 1/R compensates for the suppression of Planck scale. It is notable that such a mass 
scale is not observable by the current experiments due to very limited resolutions. In order to 
avoid of any divergency in the production cross sections arising from summation over the infinite 
number of KK modes, an ultraviolet cutoff scale A uv is essential to regulate the processes. The 
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extra dimension model is a low energy effective theory which is valid below the onset of quantum 
gravity scale, denoted by the scale Mg. Throughout this analysis, the cutoff scale of the effective 
theory A E y is chosen to be equal to Mg■ In general, Mg is different from the Planck scale in 
the presence of extra dimensions M E but it is related to Mjj according to the following relation 

mm 


M S = r(-^) 


Hed ,1 1 /( n sn + 2 ) 


x Me, 


( 2 ) 


The cross sections of processes in the ADD model are usually parametrized using the parameter 
rjc which is equal to F/Mg where F is a dimensionless parameter which takes different forms in 
different conventions including the GRW [23], HLZ [23] and Hewett [25] conventions. In the HLZ 
conventions, F is expressed as a function of Mg and number of extra dimensions: 


F = 


tog(-f-) if n ED = 2 

if n E D > 2 , 


(3) 


where y/§ the center-of-mass energy of the hard process which is approximately equal to the di- 
Higgs invariant mass in this analysis. In the GRW convention, F is equal to one and the scattering 
amplitude of graviton mediated processes can be parametrized in terms of a single parameter At 

m- 


A = ~^T , for n ED > 2 , (4) 

where T is a function of energy-momentum tensor. The GRW and HLZ conventions can be related 
using \ll = F x Aj.. In this work, the results are presented in both conventions. 

There are different possibilities which can be used to study the effects of large extra dimensions. 
Some experimental tests of the ADD model are mentioned here. 

Gravitational law: the Newton’s gravitational force will be modified in the ADD model framework 
at distances shorter than the size of the extra dimensions. At 95% CL, the size of the extra 
dimension above 37^m has been excluded. This is corresponding to an exclusion of M E below 1.4 
TeV for two extra dimensions [28], 29]. 

Collider experiments: as mentioned before, large extra dimension leads to direct production of 
gravitons at particles colliders as well as enhancements in the cross sections of some SM processes 
due to virtual graviton exchange. Experimental limits on the extra dimensions have been set by 
different experiments including the HERA [30] (31], LEP 32. [33l 134. 35j and Tevatron [36] [57]. 
At the LHC, the ADD model has been probed in diphoton, dilepton, monophoton, and monojet 
channels in both CMS and ATLAS experiments [351 ESI 501 ?, E3 (3L 321 S3. 33]- The most 
stringent collider limits on A E come from the LHC run at the center-of-mass energy of 8 TeV 
from dilepton and monojet events which are 4.0 TeV and 3.74 TeV, respectively [42, SSI IS] • 
Another interesting signature of the large extra dimensions at collider experiments is the black 
hole production [45, 06]. 

Cosmological and astrophysical constraints: cosmological and astrophysical observations provide 
strong bounds on the large extra dimension model parameters. Star cooling, 7 ray diffusion 
and universe expansion during the big bang nucleosynthesis are examples of astrophysical and 
cosmological implications by which the ADD model are constrained. More details can be found 

in [33 SUES[50]. 

There are other studies on the consequences of the large extra dimensions in the electroweak 
precision test, neutrino physics etc in the literature m E3 ES E3, 55 , SO]. 
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Figure 1: SM loop diagrams for Higgs pair production via gluon fusion. 


So far, the theoretical cross section of di-Higgs production in the context of large extra dimen¬ 
sion has been calculated in the 77 , e _ e + , and pp colliders IHUSSMEDI. In this work, we perform 
a detailed search for the ADD model based on the di-Higgs production in proton-proton collisions 
at the LHC and FCC at the center-of-mass energies of 14, 33, 100 TeV. All major backgrounds 
are taken into account and the effects of a CMS-like detector is considered. The jet substruc¬ 
ture technique is utilized to capture boosted Higgs bosons and to reach a reasonable background 
rejection and efficient event reconstruction. 

2 Double Higgs production in pp collisions 

The double Higgs boson production at hadron colliders within the SM has been studied in id, 62]- 
The representative Feynman diagrams for production of two Higgs bosons at hadron colliders are 
presented in Fig[l] The di-Higgs final state proceeds through gg fusion via quark loop diagrams 
and qq annihilation. The main contribution to the total production rate comes from the loop 
diagram involving mostly the top quark in the gg fusion. Due to larger parton distribution 
functions of the gluon and very small Yukawa couplings of the Higgs boson with light quarks, the 
dominant contribution of the di-Higgs production comes from gluon-gluon fusion involving the 
triangle and box diagrams. The total cross section of di-Higgs at next-to-leading (NLO) order 
calculated assuming the top quark mass M t =173.1 GeV, bottom quark mass Mb = 5 GeV, Higgs 
boson mass Mh = 125 GeV, af°(M§) = 0.13939 and a^ LO (M^) = 0.12018 at y/s = 14 TeV, 
33 TeV and 100 TeV are 33.89 fb, 207.29 fb and 1417.83 fb, respectively [63]. For cross section 
calculation, the CTEQ 66 [B3] PDF set is used. An interesting aspect of di-Higgs production is 
the destructive interference between the box and the triangle contributions shown in Fig[l] It is 
worth mentioning that the destructive interference is not negligible so that it leads to a reduction 
of around 50% in the production rate [^j 

Within the ADD model, the double Higgs production occurs at tree level through both gg 
fusion and qq annihilation via s-channel. The representative Feynman diagrams are presented 
in Figj2j As it can be seen, in the ADD model two Higgs bosons are produced at tree level via 
virtual gravitons exchanges. The presence of the new diagrams lead to increase the total rate of 

“The amount of reduction in the total cross section due to the interference term depends on the center-of-mass 
energy of the collision. At yfs = 14 TeV, it amounts to 50%. 


4 






Figure 2: Tree level di-Higgs production diagrams with virtual graviton exchange in the ADD 
model. 



Figure 3: Cross section of pair Higgs production at three different center-of-mass energies of 14, 
33 and 100 TeV as a function of A t in the GRW convention. It can be seen that cross sections 
tend to the SM ones as A t —> oo. 

Higgs pair with respect to the SM rate. On the other hand, due to mediating spin-2 graviton one 
expects different kinematical properties between the Higgs pair from the SM and ADD model. 
These issues will be discussed more in the next sections. 

In this work, Sherpa [B5] event generator is used to generate the di-Higgs events and to calculate 
the cross sections in the ADD model. Figure [3] shows the calculated cross section of Higgs boson 
pair production at three different center-of-mass energies of 14, 33 and 100 TeV as a function 
of ADD model parameter At in the GRW convention. As expected, the total production cross 
section of two Higgs bosons grows significantly with respect to the expectation of the SM at the 
three center-of-mass energies. Due to larger phase space and PDFs, the cross section increases 
with increasing the center-of-mass energy. Because the gg and qq interactions with gravitons are 
suppressed by the Planck scale in 4 + tied dimensions, the cross section is expected to decrease 
with increasing the ADD model scale A t hence the cross section goes to the SM value when 
A t —> oo. 

Figure [4] shows the ratio of di-Higgs cross section from qq annihilation and only bb to the total 
cross section versus the ADD model parameter A t at three collision energies of 14, 33, and 100 
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Figure 4: Ratio of partonic cross section over the total cross section for the quark anti-quark 
initial state. 



Figure 5: Ratio of partonic cross section over the total cross section for the gluon-gluon initial 
state. 

TeV. The contribution of the gg fusion is shown in Fig(5} As expected the main contribution 
is coming from the gg fusion with the amount of more than 60% and 90% at A t = 2 TeV at 
the center-of-mass energies of 14 and 100 TeV, respectively. With increasing the center-of-mass 
energy, the contribution from gg fusion is increased. It is interesting to note that again with 
increasing the center-of-mass energy, the b-quark parton distribution function is increased which 
leads to larger contribution from bb annihilation at larger energies. 

We close this section by mentioning that the ADD model leads to produce considerable number 
of Higgs boson pairs in pp collisions. We will see that the increment in number of Higgs pairs in 
particular occurs at the large invariant mass of the two Higgs system Mhh ■ Such an effect will be 
used as a tool to search for the ADD model and constraining the model parameters at different 
collision energies. The details of the analysis are described in the next sections. 

3 Analysis Method 

In this section, we explain the analysis procedure that is followed for generating and analyzing 
the signal of ADD model in pp collisions at the LHC with y/s = 14 TeV and FCC. Throughout 
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this analysis, we use a CMS-like experiment characteristics for simulating the effects of detector 
and similar statistical tools for obtaining the exclusion limits as the CMS experiment. We only 
focus on the Higgs bosons decay into bb pairs. This leads to have a final state containing four jets 
originating from the hadronization of b-quraks. 

3.1 Event generation 

The ADD signal events are generated using Sherpa version 2_1_1 [65] in the GRW convention. 
Sherpa also performs parton showering and hadronization processes. The background with most 
similarity to the signal which can be interpreted as the irreducible background comes from the 
SM di-Higgs explained previously. The SM di-Higgs event generation is done with the MadGraph 
5 [66, 67] and Pythia [68] is used to perform parton showering and hadronization. The remaining 
main background processes are QCD multijets, Zbb ZZ , ZH, tt, single top, W + W~, W + jets 
which are generated using Sherpa including the showering and hadronization. The contribution 
of multijet QCD background is difficult to reliably estimate due to large production rate. In 
reality, the determination of the QCD multijet contribution requires computing resources and/or 
employing a data-driven technique which is beyond the scope of the current work. We generate 
several QCD multijet samples in various bins of invariant mass of the final state partons with 
large amount of events in each bin. 

A full and real detector effects simulation must be performed by the experimental collaboration 
however we use Delphes [69] as the tool to estimate the response of the detector. It considers a 
modeling of the CMS detector performances as explained in m- In this study, the effects of 
pileup and underlying events are not taken into account. 

Finally, we should mention that the signal samples are generated with Sherpa in the GRW 
convention for the various values of the ADD model parameter Ay and all simulated samples are 
generated for the three possible scenarios of the center-of-mass energies of proton-proton colliders, 
14 TeV, 33 TeV and 100 TeV. 


3.2 Analysis details 

We perform the analysis of the simulated events on the stable final state particles. The selection of 
the events is designed in such a way to find the HH events with subsequent decay of HH —> bbbb. 

Before going further, one of the special characteristic of the signal events, which leads to 
employ particular strategies for events reconstruction and selection, is considered . In Fig[6j the 
normalized distribution of between two bottom quarks coming from the decay of each Higgs 
boson for two samples of signal with A t = 3,5 TeV and some backgrounds is shown. As it can 
be seen in Fig |6j signal events tend to reside at very small values of A R contrary to the SM 
backgrounds. 

In the events of signal with very large di-Higgs invariant mass, Higgs bosons are Lorentz- 
boosted particles which decay differently from the topological point of view compared to the 
Higgs bosons which are produced almost at rest. The angular separation of a bb pair coming from 
a Higgs boson can be approximated as: 


A %- 


rriH 


y/ x(l — X ) Pt 


(5) 


li Ait is the angular separation of b and b quarks in the r/ — ip plane which is defined as: A R b i = 
y/iVb - Vb ) 2 + (< fr > - M 2 
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Figure 6: Normalized distribution of the A R between two bottom quarks coming from decay of 
each Higgs for the two samples of signals and the potential backgrounds. 


where Pt is the transverse momentum of the Higgs boson, x and 1 — x are the momentum fractions 
of the b and b quarks. The larger Higgs Pt the smaller angular separation of bb pair. Figure [7] 
shows two dimensional plots of A R b i versus the Higgs boson Pt and the hardest b quark Pt in a 
Higgs boson decay for the SM process of pp —> HH —>• bbbb. The plots of Figure [7] confirm that as 
we go to the boosted region (events with large transverse momentum of Higgs or large transverse 
momentum of b quark), the A R b i decreases. It means that the boosted Higgs bosons produce one 
collimated jet with substructure. This can originate from two reasons: first is that the decaying 
Higgs boson has an energy several times larger than the Higgs boson mass in the laboratory frame 
and second reason is the difference between the mass of Higgs and b quark is large. 

As the di-Higgs invariant mass is an important quantity which will be used to separate ADD 
signal events from the backgrounds, A R b i is also presented versus the di-Higgs invariant mass 
for signal and SM di-Higgs events in Figure [8j Obviously, the ADD signal events prefer to be 
distributed in the large invariant mass region of the di-Higgs system which is not the case for 
the SM di-Higgs events. Another observation from the Figure [8] is that with increasing the di- 
Higgs mass of the signal events, A R b i decreases and falls even below 0.4. More accurately, the 
SM di-Higgs events are distributed at di-Higgs invariant mass smaller than 1 TeV and peak at 
A R =1.5 while for the signal, di-Higgs invariant mass peak at values greater than 3 TeV and 
around A R = 0.3. 

The kinematics of the Higgs boson decay products is categorized by two types of event topolo¬ 
gies. The first category consists of Higgs boson pairs which are produced near the threshold. 
In this type of events (normal events), each parton is matched to a single jet. In the second 
category, Higgs bosons are produced with a large Lorentz boost, resulting in collimated jets that 
might cluster into one jet. These events are referred as boosted events which need to be treated 
differently from the normal events. 

3.3 Boosted jet reconstruction 

As it has been mentioned previously, reconstruction of b-jets specially for the signal events due to 
the existence of highly boosted Higgs bosons are crucial. According to Eqj5] if the Higgs bosons 
have transverse momentum larger than Pt > 500 GeV and if approximately both b-jets carry 
the same fractions of the Higgs boson momentum, the angular separation between two b-jets is 
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P T b [GeV] P T H [GeV] 


Figure 7: The transverse momentum of the hardest bottom quark and the Higgs boson transverse 
momentum as a function of A R for the SM di-Higgs events. 



M HH [GeV] M^GeV] 

Figure 8: Invariant mass of di-Higgs as a function of the A R between the each b-jet pairs for the 
signal with A t = 5 TeV (left) and for the SM di-Higgs events (right). 
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A R < 0.5. As a consequence, common jet reconstruction which usually is done with the cone size 
of A R = 0.4 — 0.5 would not be applicable for the most of signal events. Therefore, an alternative 
way of fat jet algorithm is used m for these boosted events. 

Now, the jet substructure analysis is described together with its application on our signal with 
four boosted b-jets in the final state. We reconstruct the fat jets using the Cambridge/Aachen 
(CA) jet algorithm [72,, |73J assuming special jet cone size of R = 1.2. Then to identify the boosted 
Higgs bosons, the procedures described in the fat jet reconstruction algorithm m is performed 
as the following. First, large radius or a fat jet J is split into two sub-jets Ji and J 2 with masses 
to./, > mj 2 . Then a significant mass drop of mj 1 < with /i = 0.667 is required, yuo is 

an arbitrary value that shows the degree of the mass drop. In order to avoid the inclusion of high 
Pt light jets, two sub-jets have to be symmetrically split by satisfying: 

min ( P T,J^ P T,J 2 ) 


rn 


< 2/cut 


( 6 ) 


where P^ ^and P^ j 2 are the square of the transverse momentum of each sub-jet and y cu t is 
one of the parameters of the algorithms which determines the limit of asymmetry between two 
sub-jets. Finally, if the criteria in the above steps are not fulfilled, we take J = J\ and return to 
the first step to perform decomposition. All the above steps are followed by a filtering in which 
a re-clustering is performed with the radius of Rjut = min( 0.3, ARj lt j 2 /2) which selects at most 
three hardest jets. This is a useful step to remove the contributions from pileup and underlying 
events m- In the analysis, the two hardest objects are required to be tagged as b-jets while the 
third one can be possible radiation of the two b quarks. 

It has been shown the best values for the algorithm parameter are y md = 0.67 and y cu t = 0.09 
m and the best performance for clustering when it deals with jet substructure is C/A algorithm. 
It is worth mentioning here that tighter value of Hmd would not be useful significantly [75]. The 
algorithm explained above for reconstruction of the boosted objects has been implemented in the 
FastJet3.1.1 package m using that we perform the analysis to find the two Higgs bosons in the 
final state. 

In this analysis, first the jets coming from signal or backgrounds are reconstructed with the 
anti -hr algorithm with the cone size of R = 0.5 then if two jets with Pp > 250 GeV are found in 
the event, the fat jet algorithm is applied. Otherwise, the event is treated as normal event. 

As for the b-tagging efficiency and mis-tag rates, similar numbers as the CMS experiment are 
used. The data driven efficiency for the b-jet identification indicates that the efficiency is as large as 
60% to 80% [77]. In our analysis, a b-tagging efficiency of 70% for jets with transverse momentum 
larger than 30 GeV and in the pseudorapidity range of |ry| < 2.5 is assumed. Mis-tagging rates of 
10 % and 1% for the c-jets and for the light jets are considered m- 


3.4 Higgs bosons reconstruction 

For reconstruction of the Higgs bosons in the final state, a y 2 algorithm is utilized to determine 
the correct assignment of b-jets to Higgs bosons candidates. It relies on the Higgs boson mass 
and other kinematics properties as constraints. All possible permutations for four or more b-jets 
are tried and the permutation with minimum x 2 is used to reconstruct both Higgs bosons. The 
X 2 which is run over the events containing at least four b-jets with Pp b > 30 GeV and |? 7 | < 2.5 is 
defined as: 

X 2 = (Mi, - M h ) 2 + (M kl - M h ) 2 + A R% + AR 2 kl (7) 

where M$j, M k i are the invariant mass of the b-jets pairs and Mh = 125 GeV. As mentioned above, 
the best combination of b-jet pairs is the one with minimum x 2 out of three possible combinations. 
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Figure 9: Normalized distributions of the reconstructed Higgs bosons for the signal and back¬ 
grounds after applying the fat jet algorithm and detector effects. 


Figure [9] shows the invariant mass of the b-jet pairs after applying the fat jet algorithm and a 
CMS-like detector effects. As it can be seen, the reconstructed Higgs bosons from signal have 
a very good resolution on the mass spectrum. For the sake of comparison, the reconstructed 
mass distributions of some background processes are also shown in Figj9] Because of the efficient 
performance of the fat jet algorithm which leads to a better resolution on the Higgs boson mass 
spectrum with respect to all backgrounds, imposing cut on the invariant mass of each b-jet pairs 
can suppress significant amount of the backgrounds keeping signal events at a maximum level. 

According to discussions in the first section, we expect a continuous enhancement in the rate 
of the signal events due to contribution of the very close Gkk modes of gravitons. This effect 
manifests itself mostly in the high invariant mass region of di-Higgs events where the number of 
excited modes of the Gkk are much larger. This effect has been shown previously in Fig(8j the ADD 
signal events have very large di-Higgs invariant mass while the SM di-Higgs events are distributed 
at lower invariant mass with respect to signal. Such a discriminating feature is used in the next 
section to set limits on the ADD signal model parameters. Figure [TO] depicts the invariant mass 
of the two reconstructed Higgs bosons for the signal and different sources of backgrounds. This 
plot shows the behavior of the signal and potential backgrounds in the final state mass spectrum. 
In this plot, no cut except for the acceptance cuts are applied. 

3.5 Event selection 

To select the signal events, we require events to have exactly four b-tagged jets with Pt > 30 GeV 
and \r]\ < 2.5 and no isolated lepton with Pt > 10 GeV. The missing transverse energy is required 
to be less than 10 GeV. All these requirements are denoted as cut-1. Additional cuts are applied 
for further backgrounds suppression. Due to good resolution on angular separation between b-jet 
pairs coming from each Higgs boson, we require each b-jets pair to have A R < 0.6 (cut-2). Such 
a tight cut reduces the contributions of non-boosted background events. A mass window cut of 
100 GeV to 150 GeV is applied on the invariant mass of each reconstructed b-jets pair (cut-3) 
which suppresses the contribution of the backgrounds with no Higgs boson. We should mention 
here that triggering of the events with purely hadronic final states needs careful attention however 
requiring four energetic jets in the event has good enough efficiency at the LHC experiments. 

Finally, we present the event yield at the center-of-mass energy of 14 TeV with 300 fb -1 of 
the integrated luminosity in Table [l] In this table, the number of remaining events are presented 
in a signal region of Mhh >1-4 TeV for signal and main backgrounds. Table [l] explicitly shows 


11 




















Figure 10: Invariant mass of the Higgs pair for the signal sample with the parameter of At = 5 TeV 
and other sources of the backgrounds, distribution is normalized to the number of the each samples 
separately. 



ADD,A t = 6 TeV 

SM,HH 

QCD 

Z bb 

zz+ww 

ZH 

tt 

Cut on Mhh 

M hh >1-4 TeV 

cut-1 

3.86 

13.46 

5.6e+5 

1.2e+5 

26.47 

4.34 

40.70 

cut-2 

3.51 

1.96 

13.48 

31.0 

0.01 

0.22 

0.00 

cut-3 

2.28 

0.55 

0.02 

0.30 

0.00 

0.00 

0.00 


Table 1: Number of survived events of ADD signal and different backgrounds after applying sets 
of cut-1, cut-2, and cut-3 described in the text in the signal region of Mhh > 1-4 TeV for the 
yfs = 14 TeV and with 300 fb~ l of integrated luminosity of expected data from the LHC. 


that in the high invariant mass of di-Higgs, i.e. the signal region, the contribution of the reducible 
backgrounds are almost negligible and the only main source of the background comes from the SM 
Higgs pair production. As it can be seen and shown in m, the combination of jet substructure 
techniques, b-tagging requirement, and invariant mass cuts reduces the contribution of the QCD 
multijet background negligible. 

The cut on the di-Higgs invariant mass needs to be optimized to achieve best limits on the 
ADD model parameter which will be explained in the next section. 

4 Limit calculation 

In this section we present the statistical procedure that we use to obtain the expected limits on 
the ADD model parameter. As we mentioned formerly, the invariant mass of the Higgs pair is 
an effective observable that we use to set limit on the signal cross section and then translate the 
limit on the model parameters in the absence of any indication of the ADD model signal. A single 
bin counting experiment in the signal dominant region (high invariant mass region) is used to set 
the limits. We begin with a Poisson distribution as the probability of measuring N events in the 
signal region: 
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P(N\ a add eC,B) = e ~^+cr ADD eC) ( B + °add£C) n (8 ) 

where ctadd, £, £ and B are signal cross section, signal efficiency, integrated luminosity and 
expected number of backgrounds. In the above equation, a add is taken as a free parameter to 
enable us to consider different ADD signal production cross sections. To obtain the number of 
expected background events B and the signal efficiency e, we rely on the Monte Carlo simulations. 
At confidence level of 95%, an upper limit on the signal rate ctadd is obtained by integrating over 
the posterior probability as the following: 


t 95 % 


0.95 = 


fo P{N\ctaDD£ £,B)d(JADD 


(9) 


/ 0 °° P(N\ a add £ B)d(JADD 

In order to extract the expected limit on the ADD signal cross section, one has to solve the Eqj9] 
for <j 95 % under the assumption of N = B after inserting the proper inputs for the background 
expectation, signal efficiency, and the integrated luminosity. 

In the first step of limit setting, we have to determine the signal dominant region. Therefore, 
the di-Higgs invariant mass cut which determines this region is optimized in such a way that gives 
the best exclusion limits on the model parameter At- This can be reached by minimizing the 95% 
CL expected limit on the signal cross section. Figure [TT] shows the calculated expected limit at 
95% CL on At as a function of invariant mass cut of di-Higgs. As it can be seen in Figjll[ with 
increasing the cut on di-Higgs mass the exclusion limit on A^ is maximized at the cut on di-Higgs 
mass of around 1.4 TeV with an integrated luminosity of 300 fb -1 . Thus, we take the mass cut 
of 1.4 TeV as the optimized value to introduce the signal region. It has to be mentioned that in 
the optimization process no systematic uncertainty is included. It is worth mentioning that the 
optimized cut on the di-Higgs invariant mass varies with the integrated luminosity. 

We calculated the signal efficiency after applying the set of cuts described in the previous 
section. This efficiency has almost a flat behavior against the model parameter A t ■ The mean 
value of the signal efficiency is taken to calculate the limit. It is found to be equal to 17%. The 
uncertainty on the efficiency originating from statistical fluctuations and a 5% uncertainty due to 
the fluctuation of efficiency for different A t are considered. An overall uncertainty of 30% on the 
number of background events in addition to the statistical uncertainty is considered. 

Figure 12 shows the expected limit at 95% CL as a function of model parameter A t including 
the uncertainty bands. The theoretical cross section of ADD signal also is presented for compar¬ 
ison. The 95% CL expected upper limit on signal cross section in the signal region is found to 
be 0.09 fb for an integrated luminosity of 300 fb -1 of data. It leads to an expected lower limit 
on A t to be 5.1 TeV. These bounds are in a reasonable agreement with the 3 a exclusion limits 
calculated in m where no object reconstruction and identification and detector effects have been 
considered. 

Similar analyses are performed for higher center-of-mass energies of future planed proton- 
proton colliders with the integrated luminosities of 300 fb -1 and 3 ab” 1 . The results for expected 
limits on the A t in GRW convention and on ( tied , Ms) in the HLZ convention are summarized 
in Tables [2] and [3j respectively. Moving to larger center-of-mass energy of the collisions leads to 
increase the lower limit on At- The limit is extended up to around 24 TeV at the collision energy 
of 100 TeV. Using more integrated luminosity of data would lead to improve the limit on the 
model parameter. Increasing the integrated luminosity by a factor of 10 changes the lower limit 
on A t from 5.1 TeV to 6.8 TeV at y/s = 14 TeV. 

At the end of this section, it must be mentioned that including the other decay channels of 
the Higgs bosons ( 77 , t + t~ , ZZ , WW) would improve the limits considerably. 
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Figure 11: The 95% CL expected limit on At as a function of cut on the invariant mass of two 
reconstructed Higgs bosons. 
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Figure 12: The 95% CL expected limit on the ADD signal cross section as a function of A t 
including the uncertainty bands. The theoretical cross section of ADD model is also presented for 
comparison. 


IL = 300 fb~ l 
center-of-mass energy 

A t [TeV](GRW) 

M s [.TeV](HLZ) 
n=3n=4n=5n=6n=7 

14 TeV 

5.1 

6.1 5.1 4.6 4.3 4.1 

33 TeV 

10.5 

12.5 10.5 9.5 8.8 8.3 

100 TeV 

23.6 

28.1 23.6 21.3 19.8 18.8 


Table 2: 95% CL expected limit on the parameters of the model in both GRW and HLZ conventions 
for the 300 fb of integrated luminosity of data. 
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IL = 3 ab~ l 
center-of-mass energy 

A t [TeV](GRW) 

M s [TeV]{HLZ ) 
n=3n=4ra=5n=6n=7 

14 TeV 

6.8 

8.1 6.8 6.1 5.7 5.4 

33 TeV 

13.4 

16.0 13.4 12.1 11.3 10.7 

100 TeV 

28.7 

34.1 28.7 25.9 24.1 22.8 


Table 3: 95% CL expected limit on the parameters of the ADD model in both GRW and HLZ 
conventions for the 3 ab~ 1 data. 


5 Di-Higgs angular distribution 


An interesting feature of di-Higgs production from the large extra dimensions is the quite different 
behavior of the angular distribution of the final state with respect to the SM backgrounds. In 
general, final state particles coming from the exchange of gravitons with spin 2 should have 
different shape from the final state particles from the exchanges of photon, Z-boson or Higgs 
boson. On the other hand, using fat jet algorithm enabled us to have very good resolution on 
angular separation. Therefore, angular distribution of the Higgs boson pairs can be used as a 
powerful observable to distinct the ADD signal from the SM backgrounds to set limits on the 
model parameters. 

The shape of the angular distribution of di-Higgs, which is an interesting feature of the ADD 
model, is shown in Fig|13[ The angular distribution of the SM di-Higgs is presented for comparison. 
In this plot, 0{H\,H2) is the angle between the directions of the momenta of the final state 
Higgs bosons. The distribution of cos 0(H\, H 2 ) is plotted for ADD di-Higgs production from qq 
annihilation and gg fusion separately. As it can be seen, the angular distribution of the signal 
events from gg fusion has quite different behavior from qq. The two Higgs of ADD events produced 
from gg prefer to fly mostly perpendicular to each other while the two Higgs bosons in the ADD 
events from qq fusion tend to be produced at the angles of approximately ±7 t/4 . Detailed analytical 
explanations of the angular distributions of di-Higgs from the ADD model in e + e _ collisions and 
77 collisions can be found in m and [79]. Similar explanations are valid for the hadron colliders 
with the initial states of qq and gg are expected to be like e + e _ and 77, respectively. According 
to Figjl3], the SM di-Higgs distribution is quite flat and have different shape from the ADD signal 
events. It is worth mentioning here that as discussed in section 1, due to larger gluon PDF at 
larger center-of-mass energies the contribution of gg fusion in ADD signal production is increased. 
It has been shown in Fig. [5j 

The CMS and ATLAS collaborations [80j E] have used a variable of X = exp(\yi — 2 / 2 1) to 
search for the contact interactions in di-jet events. The rapidities of the two jets are denoted by 
2 / 15 2 / 2 - The rapidity y is defined as log(^j^) where E is the energy and p z is the z-component of 
the momentum of a given particle. The advantage of the rapidity difference is that it is a boost 
invariant quantity. Now, we use the X = exp(\yu 1 — yH 2 \) distribution to probe the effects of 
ADD model instead of cos 0(H\, H 2 ). 


Figure 14 shows of X = exp^yHi — Ph 2 1) f° r main SM background di-Higgs and the ADD 
signal with A t = 3 TeV. Considering only the SM di-Higgs as the main background, we set a 
limit on the ADD parameter using this angular distribution. We perform the analysis with an 
integrated luminosity of 300 fb -1 . We define a y 2 over the X distribution as: 


x 2 = 


£ 

i =0 


{■^sig+bck j 




( 10 ) 


syst 
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Figure 13: The distribution of cos 0(Hi, H 2 ) for the SM di-Higgs background and the extra di¬ 
mension signal at parton level. 



Figure 14: The normalized distribution of exp (jjh^ ~ Vh^) f° r the ADD signal events and SM 
di-Higgs production at parton level. 


where the numerator shows the square number of signals and A s t a t and A sys t represent the sta¬ 
tistical and systematical uncertainties. To calculate the x 2 , we use the same event selections as 
before and we consider the events in the signal region which is determined in the previous section. 
Figure [l5] shows the \' 2 /n as a function of At, where n indicates the number of degree of freedom. 
The dashed-line shows the value of x 2 / n which corresponds to the 95% CL. The limit on A t 
using this observable is found to be 6.98 TeV which is higher than the value that we obtained 
using the invariant mass. To check the effect of systematic uncertainties, we considered 20% and 
40% overall systematic uncertainties which are shown as dashed lines. Including 40% systematic 
uncertainty leads to loosen the limit on A t around 200 GeV. A significant improvement could be 
achieved using the distribution of the rapidity gap of two-Higgs bosons which amounts to around 
1.8 TeV with respect to the mass spectrum analysis of di-Higgs. 

This leads to conclude that the di-Higgs final state would be a promising channel to search 
for the large extra dimension effects at the hadron colliders. In particular, the usage of angular 
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Figure 15: The y 2 /n as a function of At, where n represents the number of degrees of freedom. 
The dashed line shows the value of x 2 /n which corresponds to the 95% CL. 


distribution would lead to stringent bounds on the model parameters. 

6 Summary and conclusions 

Double Higgs boson production at hadron colliders provides the possibility to probe not only the 
Higgs self-coupling and Higgs couplings with the SM particles but also it enables us to search 
for the effects of new physics beyond the SM. In this paper, the double Higgs production at the 
LHC and future circular collider (FCC) with the center-of-mass energies of 14 TeV, 33 TeV, and 
100 TeV is used to search for the effects of large extra dimensions. The analysis is only based 
on the most probable final state i.e. pp —> HH —> bbbb which is of course a challenging channel 
due to large QCD background and triggering the events. The tail of the invariant mass of the 
two Higgs bosons is affected due to the virtual gravitons exchange. In addition to the di-Higgs 
invariant mass, the angular distributions of the final state Higgs bosons (and consequently the 
decay products) have a quite different behavior with respect to the SM irreducible background due 
to the exchange of spin 2 gravitons. We perform a comprehensive Monte-Carlo simulation analysis 
taking into account the main backgrounds and consider a CMS-like detector effects using the 
Delphes package. To reconstruct the signal candidate events efficiently and reasonable background 
rejection, jet substructure techniques are employed to capture the signal events which are boosted 
objects in the final state. Then we obtain the expected limits on the model parameters using the 
invariant mass and the angular properties of the final state particles. Depending on the number 
of extra dimensions, the effective Planck scale is limited up to 6.1, 12.5, 28.1 TeV at the proton- 
proton collisions with the center-of-mass energies of 14, 33, and 100 TeV, respectively. Further 
improvement of the analysis is possible by including other decay modes of the Higgs bosons such 
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as 77, WW, ZZ, t + t~ . 
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